Determination of Mueller matrix of an optical element with Simon-Mukunda polarization gadget 
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We have described a novel way to determine the Mueller matrix of any optical element by using projection 
method. For this purpose, we have used two Simon-Mukunda gadgets to obtain projection matrix directly from 
the experiment. Mueller matrix has been determined using the experimentally obtained projection matrix for 
three known optical elements namely free space, half wave plate and quarter wave plate. Experimental matrices 
are in good agreement with the corresponding theoretical matrices. 



I. INTRODUCTION 

To characterise the polarization changes in a given light 
beam due to any optical element, one needs to find its Mueller 
matrix. From the Mueller matrix, we can obtain the informa- 
tion about the retardation introduced, the diattenuation and 
the depolarization of a light beam after passing through any 
optical element [ 1 1 2 1. Determination of Mueller matrix have 
diverse applications in various fields such as tomography 131 . 
micro-electronics |4i|, geo-science [5|, astronomy J5), study 
of polarization changes in vortex beams [7|, liquid crystals 
and characterisation of undesired modulation introduced by a 
spatial light modulator |S). Along with these, determination 
of Mueller matrix in optical coherence tomography gives the 
birefringence and the orientation of optic axis with respect to 
the incident beam for a given sample 191 . 

Simon-Mukunda (SM) polarization gadget, an universal 
gadget, can be used to generate any polarization state on the 
Poincare sphere with minimum number of optical elements 
IfTOlfTTl . This gadget consists of one half wave plate (H) and 
two quarter wave plates (Q); the polarization states can be ob- 
tained by changing the rotation of their fast axes appropriately 
IfTSl . One can use SM gadget in different configurations such 
as Q-H-Q, Q-Q-H and H-Q-Q. However, we have used Q-H-Q 
configuration only. The principle of SM gadget for determin- 
ing the polarization states has been used in variety of applica- 
tions in recent years |fT3TfT5l . In the depolarization studies of 
light beams due to polar meso-spheric cloud, Lidar can intro- 
duce additional polarization changes to the receiving signal. 
To compensate these changes, Hayman et.al. have used the 
principle of SM gadget |fT3l . Schilling et. al. presented a the- 
oretical scheme to determine the higher order correlations in 
a single light beam by introducing proper unitary transforma- 
tions to the input state with SM gadget [14]. While measur- 
ing the Pancharatnam phase of a polarization state with tech- 
niques of interferometry and polarimetry, Loredo et. al. have 
used SM gadget [15]. However, in the present study, we have 
used it to find out the Mueller matrix of a given optical ele- 
ment. We have used two SM gadgets; one for the generation 
of input polarization states and another for obtaining the pro- 
jections of each output state on a given set of four input states. 
To the best of our knowledge, we are the first to demonstrate 
the use of this gadget to determine the Mueller matrix. 



In conventional approach, Mueller matrix is determined 
by measuring the output Stoke's parameters of a given light 
beam, corresponding to four different input Stoke's parame- 
ters. These parameters are found by measuring the intensities 
with different polarizations (horizontal, vertical, left and right 
circularly polarized and both diagonal). This method is also 
known as successive probing method [ 16] and it needs 24 in- 
tensity measurements with known input Stoke's parameters. 
In Mueller matrix imaging polarimetry (MMIP), one can ob- 
tain Mueller matrix by taking 16, 36 or 49 images ifTTl . In 
one of the earlier studies, to obtain the Mueller matrix, Az- 
zam presented a technique called photo polarimetry that uses 
two quarter wave plates (one in polarizing optics and another 
in analysing optics) rotating at speeds to and 5oj fl8l . Alterna- 
tively, the present study provides a novel and simple method 
to determine the Mueller matrix of an optical element using 
two SM gadgets. The theory of this method is discussed in 
section 2 and experimentally demonstrated in section 3. 



II. THEORETICAL ANALYSIS 

The Stoke's vector describes the complete polarization state 
of an electromagnetic wave. When an optical beam with 
Stoke's vector (S) passes through an optical element, the out- 
put Stoke's vector (S) can be written as 

S = MS (1) 

where S and S are four element column vectors and M is a 4 x 
4 matrix defined as Mueller matrix. We form two matrices Q. 
and Q by arranging the four input and output Stoke's vectors 
as four columns 

n=[S 1 S 2 S 3 S i ], Q = [SiS 2 S 3 54]. (2) 

Using Eq. ([!]), the relationship between Q and O can be writ- 
ten as 

Q = MQ.. (3) 

The four input states forming Q, must be on the surface 
of the Poincare sphere. The projection of each output state 
§ i, §2, §3, §4 on input states Si, S2, S3, S4 gives the 16 
non-negative real numbers which form the projection matrix 
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A. These elements are given by 

1 1 4 

A^ = -(S/S~=-2>JS? (4) 

where i, j = 1 to 4 and S a denotes the a th element of i th Stoke's 
vector and Ay is the projection of 5, state on S j state. There- 
fore, the projection matrix 

A = -Q r Q = -Q r MQ. (5) 

2 2 

As a consequence, Mueller matrix M can be written as 

M = 2(n r )" 1 AO" 1 . (6) 

To obtain Mueller matrix, one needs to choose four input 
states and determine the corresponding projection matrix. For 
the sake of completeness, we have considered three tetrahe- 
drons on the Poincare sphere as three sets of input states. The 
four vertices of each tetrahedron form a set of input states. We 
have selected these three tetrahedrons to optimise the errors in 
Mueller matrices fl6l . The optimal input states (\,qi,Uu v,) r 
for any polarimeter, should obey the following conditions. 

4 4 4 

&'' = X Mi = Z v '- = 0; ( ? ) 

1=1 1=1 1=1 



(0,-5^,-1.37^), (0,-3%, 2.91*2) and (0,-^,2.40*2) 
respectively where ko = cos -1 (0.68) = 47.16°. 



Similarly, the four vertices of third tetrahedron are repre- 
sented by 
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and the corresponding Euler's angles are (0, -f,*3), 
(0,-ifg, 2.30*3), (0, oo'^a) and (0, f,2.30* 3 ) respectively 
where k 3 = cos" 1 (0.58) = 54.54°. 



The sum of four vertices for these tetrahedrons on the 
Poincare sphere is origin, the center of mass and corresponds 
to the maximally mixed (completely un-polarized) state. Op- 
timal states for polarimeter should have minimum norm and 
condition number along with maximum determinant. Here we 
have considered Frobenius norms and condition numbers due 
to the direct dependence of errors in Mueller matrix on them 
ifTBTl . These norms and condition numbers of matrix S are de- 
termined by formulae 

||S = yj trace{S* S) (10) 



^ qiUi = 2^ = ^ v t qi = and 



1=1 



1=1 



1=1 



(8) 



(9) 



By solving these equations, one can get a number of optimal 
states for polarimetry. 

The Stoke's parameters (S],i — 1-4) representing the first 
tetrahedron are given by 
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and the corresponding Euler's angles are (0,0,0), 
(0,0, aci), (0,-y,*i) and (0, y,*i) respectively where 
cos" 1 (-0.333) = 109.47°. 



The four vertices of the second tetrahedron are represented 
by the following Stoke's parameters (Sf, i = 1 - 4) 
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and the corresponding Euler's angles are (0, -.((5,^2) 



Cf = norm(S) norm(S ] ) 



(11) 



The Frobenius norms, corresponding condition numbers and 
determinants for three tetrahedrons are (2.8279, 2.8280, 
2.8349), (4.4729, 4.4727, 4.4640) and (3.0755, 3.1218, 
3.0763) respectively. To generate four input polarization 
states (S i, Si, S3, S4) belonging to a particular tetrahedron, 
we have used first SM gadget. For example, for the gener- 
ation of a particular polarization state (£, 77, £), the fast axes' 
orientation of the wave plates in SM gadget are given by the 
subscripts ofO, H. O 1121. 



(12) 



where £,r],£ are Euler's angles corresponding to the polar- 
ization state. These angles represents the angle made by the 
frame of reference (£) from which observations are made, the 
phase lag between two orthogonal components (77) in a given 
polarization state and the rotation (£) of that polarization state 
from horizontal or vertical direction respectively. The angles 
of <2i, Hi, Q\ in Table (1) corresponding to S^ generate the 
i th state of k' h tetrahedron. 

To obtain the projections of output states from unknown 
optical element on a given set of input states, we have used 
another SM gadget. Let A* represents the projection of S ,■ 
state on S ■ state. S ,■ state is the output state, obtained after 

passing 5* state through the sample. To generate the required 
polarization state, we kept the angles in first SM gadget as 
k iHt and k iQ , . To project the output state on S k j state, we 
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TABLE I. Angles of rotation for wave plates in SM gadget 1 to generate the four vertices of three tetrahedrons. 



kept the angles in second SM gadget as 6jg i + |, 6^ Hi + \ and 
6 k n , + f ; it will generate the inverse of that 5*. 

(QeT 1 = Q i+ e (H e )-' = ff f+e (13) 

The sixteen projection elements corresponding to different 
orientations of wave plates form the projection matrix. We 
have done experiment for three sets of input states forming 
vertices of tetrahedron on the Poincare sphere and correspond- 
ing projection matrices are obtained. The reported Mueller 
matrix is the average of all three Mueller matrices correspond- 
ing to three sets of input states. 

III. EXPERIMENTAL DETAILS 

The experimental set-up for the generation of input polar- 
ization states and to obtain projection matrix for calculation 
of the Mueller matrix is shown in Fig. 1. A He-Ne laser hav- 
ing wavelength 632.8 nm and vertical polarization is used for 
this study. Laser beam is allowed to pass through a polariser 
Pi to increase the ratio of vertical to horizontal polarization 
of the laser and to set the reference for fast axis of all the 
polarizing elements. Vertically polarized light is then passed 
though the first SM gadget {Q\ H\ Q' y ) to generate four in- 
put states. These four states are further passed though a given 
sample whose Mueller matrix needs to be determined. The 
projections of resulting states from the sample on the four in- 
put states coming from first SM gadget with respect to initial 
polarization (set by P\) are obtained by using a second SM 
gadget (Qi Ho Q' 2 ) and an analyser (P2) whose transmission 
axis is parallel to the initial polarizer. To obtain this, the wave 
plates in second SM gadget are rotated in such a way that it 
will compensate the polarization changes occurred in the ini- 
tial state while generating them. The angles are obtained sim- 
ply by adding 90° to the angles used in first SM gadget (Table 
1). For example, to get the projections of first output state, we 



fix the first set of angles (?, -?, |) in first SM gadget and take 
measurements for four sets of angles obtained by adding 90° 
to the wave plate orientations corresponding to four polariza- 
tion states of the tetrahedron with second SM gadget. Thus, 
we repeat the experiment for four sets of angles corresponding 
to four input states to obtain 16 non-negative projection ele- 
ments as discussed in Eq. Q. The same procedure is repeated 
for other two tetrahedrons also. These elements form a 4 x 4 
projection matrix which is required to determine the Mueller 
matrix by using Eq. (|6j. The resultant Mueller matrix is the 
average of three matrices corresponding to input states of the 
three tetrahedrons. An optical multimeter with an accuracy of 
0.2% and resolution of 0.01 pW is used to measure the out- 
put power. We are using three known optical elements - free 
space, half wave plate and quarter wave plate as samples. We 
have used zero order quartz wave plates from Melles Griot for 
making SM gadgets as well as samples for measuring Mueller 
matrices. These wave plates are mounted on a motorized com- 
puter controlled cylindrical rotation stages with a resolution of 
0.04°. The quoted retardance tolerance of the wave plates is 
A/200 to /1/500. To begin with, we have aligned the fast axes 
of all the wave plates as vertical. To assure that the fast axes 
are perfectly vertical, we put an analyser with its fast axis per- 
pendicular to the initial polarizer. Each optical element is sep- 
arately aligned in such a way that we get a minimum power in 
the power meter. 



IV. RESULTS AND DISCUSSION 

We have experimentally determined the projection matrices 
for free space, half-wave plate and quarter-wave plate with 
their fast axes as vertical. The obtained projection matrices, 
represented by the symbols A^ e , where 1 = 1—3 

represents a particular tetrahedron. To obtain Mueller matri- 
ces for all three optical elements, we have used projection 
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FIG. 1. Experimental set-up for the determination of Mueller matrix using two SM-gadgets 



matrices in Eq. |6]). We have taken the average of three 
Mueller matrices corresponding to different tetrahedrons for 
a given sample. The average Mueller matrices from the ex- 
periment are denoted with symbols M^ ee , for free 
space (free), half wave plate (H) and quarter wave plate (Q) 
respectively. All the projection and Mueller matrices are nor- 
malized with their first element. 

To obtain the theoretical projection matrices, we have used 
Jone's matrix analysis. The projection elements represent the 
fraction of the input intensity of a light beam after passing 
through all optical elements. They have been obtained from 
the output Jone's vectors of the beam. The Jone's vector (J) 
of a light beam and corresponding intensity (/) is given by 



/ = \E X \ 2 + \Ey\- 



The output Jone's vector has been determined from the 
product of Jone's matrices of all optical components in the 
same order as shown in experimental set up (Fig. [TJ with in- 
put Jone's vector. 



output 



— Jp^ Jq-, Jh^ j< 



Q' 2 J sample J ' Q\ Jh { J Q\ J P { J input (14) 



We have used the following Jone's matrices for wave plates 
and polarizers Il9l . For a wave plate, with retardation <p and 
fast axis orientation 8, the Jone's matrix is given by 



Jr(4>,0) 



i<p ~ i$ ry up i(p 

e^ cos z (6>) + e~2 sin (9) (<?2 - e -y ) cos(#) sin(#) 
{e'i - e~z ) cos(#) sin(#) e f sin 2 (0) + e~ t cos 2 (6 i ) 

(15) 

and for a polarizer, with fast axis orientation 0, it is 



theoretical Mueller matrices have been denoted with symbols 
M< frle' '» M q '• ^ ne experimental projection and Mueller 
matrices are in good agreement with the theoretical matrices 
shown below. 



cos 2 (0) cos(6») sin(6») 
cos(ff)sin(6) sin 2 (6») 



(16) 



The sixteen non-negative and real projection elements have 
been obtained for 16 sets of orientations of wave plates in SM 
gadgets which form a projection matrix. The theoretical pro- 
jection matrices are denoted with symbols A^ e ,A^ , Ag 
where i = 1-3 corresponding to three tetrahedrons. The 
obtained theoretical projection matrices from output Jone's 
vectors are used in Eq. (|6]l to get the Mueller matrix. The 
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physically realizable or not with calculations of N-matrix dis- 
cussed in 1 20 1 . The N-matrices for three experimental Mueller 
We have verified if the experimental Mueller matrices are matrices are given by 
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The eigenvalues of the N-matrices are (-0.0117, 0.0036, 
0.0176, 2.0024), (-0.0189, 0.0039, 0.0111, 2.0055) and 
(-0.0079, -0.0033, 0.0136, 2.0071) for free space, half wave 
plate and quarter wave plate respectively. From these values, 
it is clear that all eigenvalues are real and only one eigenvalue 
is non-zero. It shows that N-matrices are hermitian and the 
Mueller matrices represent the real physical systems. Apart 
from this, we have also verified the following inequalities 



N 2 = tr(N).N, 



tr{N) = 2m Q0 and tr(MM T ) < 4m, 



no- 



(17) 



(18) 



Eq. (17 1 proves that the chosen optical elements are non- 
depolarizing. The last equation is a constrain on degree of 
polarization for a physical system, that is also satisfied in our 
case ll2"D . We have calculated the retardance for QWP and 
HWP by decomposing their Mueller matrices [ 1 ]. The retar- 
dations of HWP and QWP are 0.997tt, 0.501tt respectively. 
These errors are with in the quoted tolerances by the manu- 
facturers. 



V. ERROR ANALYSIS 

The main sources of errors in this set up are the deviations 
in fast axes of wave plates from vertical axis and errors in their 
retardance. Along with these, fluctuations in the intensity of 
light and the measurement by the power meter also contribute 
to errors in Mueller matrix of the element. We have neglected 
the error due to tilt of the wave plates which has been taken 
care in the alignment of the set up. The error analysis of the 



total system is done by the procedure followed in Il22f[24l . 
The whole experimental set up can be realized with the Jone's 
matrices and written as fl9l 



A 



(19) 



,, - / [ Jp{0p2) . Jr^RQ, Rb j). 

Jlt(<pRH, 0R5j) ■ Jr(0RQ, &R4j) ■ O . 

Jr(<Prq, 8m) ■ Jr(4>rh, Oru) ■ 
Jr{(Prq,6 ru ).J p {6 px ).L\ 



where / denotes the intensity recorded by the detector, Jp{6) 
denotes the polarizer rotated by an angle 9, Jr(4>, 0) denotes 
the retarders having retardation <p rotated by an angle 9, RH 
and RQ denote the retardation offered by HWP and QWP re- 
spectively, PI and P2 are the angles of rotation of polarizers, 
O is the object for which the Mueller matrix is determined and 
L denotes the laser. A denotes the power shown by the power- 
meter. In Jones vector study, L representing a vertically po- 
larized light with a 2 x 1 column vector and all other optics 
are 2x2 matrices. The full quantity inside the function / is 
a 2 X 1 matrix. The power meter values provide the elements 
of the projection matrix. From the different manufactures of 
the optics and rotation stages, we obtained the tolerances or 
mean error introduced by the individual elements. Consider- 
ing these errors to be uncorrected, we can write the variance 
introduced in A due to all the optical elements is | 



OA,, = 



z 



df 

dx 



(20) 
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where x's are the variables from Eq. ( 19 1. Tolerances in opti- 
cal elements and their rotation angles are 

o-Opi = 0.03° 

O-0P2 = 0.03° 
crcf> RH = 1.26° 
o-<p RQ = 1.26° 

crQ m = 0.04°, for i = 1 to 6. (21) 



Using eqs. 19pT we can evaluate the errors corresponding to 
all the elements of A matrix. After obtaining the error in pro- 
jection elements, we have introduced the 1% intensity varia- 
tion due to laser and 0.2% due to the detector. Now we do 
the same error propagation analysis with 16 variables present 
in projection matrix to obtain the error in Mueller matrix ele- 
ments as 



0~m = 



0.007 0.012 0.012 0.013 

0.012 0.021 0.021 0.022 

0.012 0.021 0.019 0.020 

0.013 0.022 0.020 0.026 



(22) 



Experimentally observed errors are close to the theoretically 
calculated errors in Mueller matrix. In our experiment, the 
maximum error obtained in individual elements while calcu- 
lating Mueller matrices are ±0.020, +0.016, +0.013 for free 
space, half wave plate and quarter wave plate respectively. 



Now we compare our results with other methods based 
on polarimetry. The maximum error in the Mueller matrix 
determined by Dev et al.[26| with simple polarimetry is 
±0.0985. Goldstein et al. l25l experimentally determined 
Mueller matrix with the maximum error of ±0.034 by using 
dual rotating retarder method. Bueno l27l obtained Mueller 
matrix with the maximum error of ±0.014 by using the liquid 
crystal variable retarders whose retardance is changed by 
applied voltage. Baba et al. determined Mueller matrix by 
imaging polarimetry that uses variable rotators and retarders 
E51 . They have obtained the maximum error ±0.0353 by 
taking 16 images and ±0.0138 while taking 36 images. We 
have taken the Mueller matrix of free space for comparing the 
error with other methods. The present method uses simple 
wave plates and takes just 16 measurements. 



These errors are within the limits as calculated in eq. 22 



VI. CONCLUSIONS 

We have discussed a new and simple technique to deter- 
mine the Mueller matrix of an optical element by using two 
SM gadgets and it's efficacy has been verified for three known 
optical elements i.e. free space, half wave plate and quarter 
wave plate. One may use this gadget in quantum process to- 
mography [30 1 and polarization state tomography |29| which 
finds diverse applications in quantum information. 
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